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A method is proposed for constructing an exact ground-state wave function of a two-dimensional 
^\ ' model with spin 1/2. The basis of the method is to represent the wave function by a product of 

0^ . fourth-rank spinors associated with the sites of a lattice and the metric spinors corresponding to 

0^ ' bonds between nearest neighbor sites. The function so constructed is an exact wave function of 

a 14-parameter model. The special case of this model depending on one parameter is analyzed in 
^ , detail. The ground state is always a nondegenerate singlet, and the spin correlation functions decay 

Q ' exponentially with distance. The method can be generalized for models with spin 1 /2 to other types 

of lattices. 
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1. INTRODUCTION 



. There has been growing interest lately in quantum spin systems with frustrated interactions. Of special impor- 
tance are models in this category for which it is possible to copstruct an exact ground state. The first example of such 
a representation is the well-known Majumdar-Ghosh model.til It comprises a chain of spins 1 /2 with antiferromagnetic 
interactions Ji and J2 of nearest neighbor and next-nearest neighbor spins, where J2 — Ji/2. The ground state of this 
"j^ ' model is two-fold degenerate and consists of dimerized singlets; moreover, there is a gap in the spectrum of excitations. 
Another example of an exactly solvable model is the one-dimensional model with bilinear and biquadratic interactions 
I and spin 1, investigated by Affleck, Kennedy, Lieb, and Tasak£3 (AKLT model). Its ground state has a structure of 
■ O the type where each neighboring pair of spins has valence bonds. It is not degenerate, the spin correlation functions 
C in the ground state decrease exponentially with distance, and there is a gap in the spectrum of excitations. This 
^ model therefore has properties predicted by HaldanetS for the one-dimensional Heisenberg antiferromagnetic model 
I I with spin 1. The valence-bond ground state is also exact for systems with many dimensions, but with spin d/2 (d is 
the coordination number of the lattice). L3 The one-dimensional AKLT model has subsequently been generalized and 
investigated in a number of papers,E3 where it has also been shown that the wave function of the ground state can be 
represented by the trace of the product of matrices describing the spin states of sites of a chain (the "matrix" form). 
These two examples are characterized by the fact that the total Hamiltonian of the model is written as a sum of cell 
' Hamiltonians (which are not mutually commuting), and the exact ground-state wave function of the total system is 
■ the eigenfunction haviiifi; the lowest energy of each cell Hamiltonian. 
T-H ' We have previouslyliS investigated an exactly solvable, one-dimensional, frustrated model with spin 1/2, whose 
0^ , properties by and large are similar to those of the AKLT model. The ground-state wave function has a special 
recursion formula, and we have shown that it can be reduced to matrix form. It must be noted, however, that both 
the recursive form and the matrix form are essentially one-dimensional constructions and cannot be extended directly 
to higher-dimensional systems. We cite Ref. |l^ in this regard, where a method has been proposed for constructing 
I an exact wave function of the ground state for models with spin 3/2 on a hexagonal lattice. The same method is 
I applicable to other systems with spin d/2. 
r-| ' In this paper we consider a class of models with spin 1/2 for which the exact wave function of the ground state can 
O , be represented in an alterpative form. In the one-dimensional case this wave function reduces to a wave function that 
O ' we have found previouslyjij but it admits generalization to higher-dimensional systems. The present study is devoted 
^ ] primarily to an analysis of the two-dimensional model. 
" '~j ■ The article is organized as follows. In Sec. we discuss the method of construction of the exact wave function for a 
rS [ one-dimensional model with s — 1/2. In Sec.]^ we formulate an exactly solvable two-dimensional model. In Sec. ^ we 
investigate the properties of this model with the aid of numerical calculations. In Sec. ^ we discuss the possibility of 
generalizing our treatment to other types of lattices. The Appendix gives a proof of the nondegeneracy of the ground 
state of the two-dimensional model in the presence of cyclic boundary conditions. 



2. ONE-DIMENSIONAL MODEL 



We have previo uslyE investigated a one-dimensional, one-parameter model containing two spins 1/2 in the unit cell 
and admitting exchange interactions between nearest neighbor spins and spins separated by two sites of the lattice. 
The exact ground-state wave function of the cyclic chain can be written in the form 
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vl-o = Tr [Ail, 2) A(3, 4) . . . A(iV - 1, N)] , (1) 

where A{2i — 1, 2 j) is a 2 x 2 matrix associated with the ith unit cell. 

Below we write the wave function hi a form more suitable for subsequent generalization to other types of lattices 
and give the general form of the Hamiltonian for which is an exact wave function of the ground state. 

We consider a chain of = 2M spins 1/2. The wave function of this system is described by the Mh-rank spinor 

* = ^f-^'^"-^, (2) 

where the indices A, /i, z^, . . . , t = 1, 2 correspond to different projections of the spin 1/2. 

We partition the system into pairs of nearest neighbor spins. The wave function can then be written as the product 
of M second-rank spinors 

^' = ^'^A'(i)^t'P(2)...^'^^(M). (3) 

We now form a scalar from Eq. simplifying the latter with respect to index pairs: 

= «'\(1)*^(2)...*"a(M). (4) 

Here subscripts correspond to the covariant components of the spinor, which are related to the contravariant compo- 
nents (superscripts) through the metric spinor 

3A,=/^= J). (5) 

*A = .gAM^-^, = .9^^*p. (6) 

The scalar function can thus be written in the form 

= ^^''{l)g^,.^-P{2)gp^ . . . ^''^M)grx. (7) 

The scalar function 'i's does not depend on the angle of rotation of the coordinate system and, hence, corresponds to 
the singlet state. 

The second-rank spinor describing the pair of spins 1/2 can be written in the form 

^'^A' ^ ct^'f'' +c,*^^, (8) 

where 5*^^ and ^'^^ are symmetric and antisymmetric second-rank spinors, respectively, and Ct and Cs are arbitrary 
constants. We know that the symmetric second-rank spinor describes a system with spin 1, so that the pair of 
spins 1/2 in this case forms a triplet. If "if^^ is an antisymmetric second-rank spinor reducible to a scalar multiplied 
by gxfi, the spin pair exists in the singlet state. Consequently, the ratio of the constants ct and Cg determines the 
relative weights of the triplet and singlet components on the pair of spins s = 1/2 and is a parameter of the model. 
In particular, for Cs — the wave function (||) contains only a triplet component, and for cj = it contains only a 
singlet component. 

In general, we can make the ratio of the constants Cg/ct different in different pairs, but to preserve translational 
symmetry, we confine the discussion to the case in which this ratio is the same in every pair. 

We note that the wave function (Q) has the matrix form (0), the matrices A{2i — 1, 2i) representing a mixed 
second-rank tensor: 



AA.(l,2) = «'\(l) = ct 



^ T:i^if^2 + I3ia2) P1P2 \ 1 

^ 1 



-aia2 --j^{oLi[32 + Pia2) 



-iaiP2-Pia2)I, (9) 



where and f3i denote the up and down projections of the spin s^, respectively, and / is the unit matrix. 

We now choose a Hamiltonian if for which the wave function is an exact ground-state wave function. To do so, 
we consider the part of the system (cell) consisting of two nearest neighbor spin pairs. In the wave function (0) the 
factor corresponding to the two spin pairs is a second-rank spinor: 

^^^{i)g^.'f''p{t + l). (10) 
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In the general case, therefore, only two of the six multiplets forming two pairs of spin 1/2 — one singlet and one 
triplet — are present in the wave function (|o|). Inasmuch as four spins 1/2 form two singlets and three triplets, the 
specific form of the singlet and triplet components present in the wave function (|l^) depends on the ratio Cs/ct- The 
cell Hamiltonian acting in the spin space of nearest neighbor spin pairs can be written as the sum of the projectors 
onto the four missing multiplets with arbitrary positive coefficients Ai, A2, A3, A4: 

i?.,+i=^A,/^'^+\ (11) 
fc=i 

where P]^^'^^ is the projector onto the missing multiplets in the corresponding cell Hamiltonian. 

The wave function (|^) is now an exact wave function of the ground state of the cell Hamiltonian with zero 

energy, because 

i?M+il*s)=0. (12) 

and Ai, A2, A3, A4 are the excitation energies of the corresponding multiplets. 

The total Hamiltonian of the entire system can be written as the sum of mutually noncommuting cell Hamiltonians: 

N 

H = Y,H,^,+u (13) 

i=l 

and since each term in (^|) yields zero in its action on ^I'^, we have 

i/|*,)=0. (14) 

The nondegeneracy of the ground state of this Hamiltonian has been rigorously proved.lll 

Since the specific form of the existing and missing multiplets in the wave function (|^) on every two nearest neighbor 
spin pairs depends on the model parameter Cs/cj, the projectors in (|l^) also depend on Cg/ct. Each projector can be 
written in the form 

-Pfc'^ = '^12^ (Sl • S2 + S3 • S4) + j}3^(si • S3 + S2 ■ 84) + j}4^Si • S4 
+ J23^S2 • S3 + j[''\si ■ S2)(S3 ■ S4) + J;^''^ (si • S3)(S2 ■ S4) 

+ (81-84) (S2- S3) + CW, (15) 

and this representation is unique for a fixed value of the parameter Cs/ct- 

Substituting the above expressions for the projectors into Eq. (|Tl|), we obtain the general form of the cell Hamilto- 
nians Hi^i^i. Inasmuch as the Hamiltonians Hi^i+i have exactly the same form for any % it suffices here to give the 
expression for Hi,2'- 

Hl,2 = -^12(81 • 82 + 83 • 84) + Jl3(8i • 83 + 82 • 84) + J1481 • 84 + J2382 • 83 

+ Ji(8i • 82)(83 • 84) + J2(8i • 83)(82 • 84) + 73(81 • 84)(82 • 83) + C, (16) 

where all volume integrals depend on the model parameter and the spectrum of excited states Ji = 
Ji{cs/ct, Ai, A2, A3, A4). In particular, for Cg — 0, choosing A2 = A3 = A4 and A1/A2 = 3, we obtain an expres- 
sion for -ffi^2 in the form 

i/i,2 = Li.L2 + i(Li.L2)2 + | (17) 

where Li = 81 + 82 and L2 = 83 + 84. 

The Hamiltonian (|l^) has the form of the AKLT Hamiltonian, a result that is not too surprising, because for 
Cs = two spins 1/2 in a pair effectively form spin 1. Note, however, that for Cs = a set of different forms of the 
Hamiltonian Hi 2 exists, corresponding to a different choice of coefficients \k- 

In general, the Hamiltonian ( |l6|) contains both bilinear and four-spin interactions. The latter can be excluded by 
setting Ji = J2 = J3 = and solving these equations for Ai, A2, A3, A4. However, since the condition Ai, A2, A3, A4 > 
0, generally speaking, is not satisfied over the entire range of the parameter Cs/ct, the simplified Hamiltonian will 
also have a ground state described by the wave function (|^) only in the region where Ai, A2, A3, A4 are positive. The 
nonzero exchange integrals J12, J13, t/14, J23 and the constant C depend only on the parameter Cs/ct- The explicit 
form of this dependence is given in Ref. in which we have also calculated the ground-state spin correlation function 
(siSj), which decays exponentially with correlation length ~ 1. 

We emphasize that the spin correlation functions {siSj) do not depend on the choice of Ai, A2, A3, A4 for a fixed 
parameter Cs/ct, because the ground-state wave function of the four-parameter set of Hamiltonians is the same. 



3 



■O o 

O <^ 



FIG. 1. Two-dimensional lattice on which the snin model is defined. 




FIG. 2. Graphical correspondence of the model wave function. The indices of the site spinors depend on the site index (not 
shown in the figure). 



3. TWO-DIMENSIONAL MODEL 



We consider an M x M-site square lattice with cyclic boundary conditions. We replace each site of the lattice 
by a square (Fig. ^ with spins s = 1/2 at its corners, making the total number of spins equal to 4M^. To avoid 
misunderstanding, however, from now on we continue to refer to these squares as sites. The wave function of the 
system is described by the product of fourth-rank spinors 



* = ]J*^"^"''"''"(n). 

n 

By analogy with (0), from Eq. (18) we form the scalar 

* = ^*'"''"'^"''"(")■9-.W5p„M„+.• 



(18) 



(19) 



where a and b are unit vectors in the x and y directions. 

The singlet wave function ( p^ ) is conveniently identified graphically with a square lattice, each site corresponding 
to a fourth-rank spinor 'i'^f^'^P (whose form is identical for all sites), and each segment linking sites corresponds to a 
metric spinor (Fig- H)- 

To completely define the wave function (p^), it is necessary to know the form of the site spinor \it-^^"^P. For this 
purpose we classify an arbitrary fourth-rank spinor, simplifying and symmetrizing it with respect to different pairs of 
indices. We have the following types of spinors as a result: 

1) a fourth-rank spinor Q^^^^p symmetric with respect to all indices; 
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S4(3) 
S2(3) 

S,(l) S^(2) 

FIG. 3. Lattice sites associated with interactions -ffi,2 and -ffi,3. 

2) three hnearly independent products of a symmetric and an antisymmetric second-rank spinor: (p^^g^p, (p^'^gf^p, 
and ip^^gpy] 

3) two hnearly independent products of two metric spinors and a scalar function: gxpgupX f-nd dx^gtipX- 
According to this classification, any fourth-rank spinor can be written in the form 

■q,xpup ^ c^gV^^p _|_ C2(^i^5i.p + 031^2" gpp + 04^3^ g^^ + c^gxpQvpXi + c65ai/5mpX2- (20) 

We note, however, that because the system of four spins 1/2 contains one quintet, three triplets, and two singlets, 
Eq. (^) still does not completely determine the form of •^^t^'^P ^ and it is necessary to determine the specific form of 
the spinors <^x^, (^2"^, and (^3'' and the scalar functions xi s-^d X2- 

Each symmetric second-rank spinor Lp^^^ describes a triplet state of the system, representing a linear combination of 
the three basis triplet functions ip^i , Vt2' ^^'^ fts- ^^ow specify nine linearly independent spinors describing 

triplet states of four spins s = 1/2: 

Au Au Au 

'Pti9iyp, Vt2 9vp, ftsgiyp, 

fti9pp, ft29pp: fa9pp, (21) 

\p Xp Xp 

Vti9pu, Vt29pyi Vtz9pv- 

The products of two metric spinors and a scalar function gx^g^pXi ^^'^ 9\v9ppX2 describe singlets states of four spins 
s = 1/2. Since there are two independent singlet functions x^i and Xs2, we have four linearly independent scalars 
describing singlet states of four spins s = 1/2: 

9\p9i^pXsii 9>.p9vpXs2, /22^ 
9Xv9ppXsi, 9\u9i^pXs2- 

As a result, the specific form of the fourth-rank spinor ^i-^^^^p [and, hence, the wave function (p^)] describing the 
system of four spins s = 1/2 is governed by 1 + 9 + 4 = 14 quantities, which are parameters of the model. 

We now choose a Hamiltonian H for which the wave function (|l^) is an exact ground-state wave function. As in 
the one-dimensional case, we seek the required Hamiltonian in the form of a sum of cell Hamiltonians acting in the 
space of two nearest neighbor spin quartets: 

-ff = ^ ^fn,n+a + ^ ^^n,n-|-b- (23) 
n n 

The first term in Eq. (|2^) is the sum of the cell Hamiltonians in the horizontal direction, and the second term is 
the same for the vertical. The cell Hamiltonians along each direction have the same form, but the "horizontal" and 
"vertical" Hamiltonians differ in general. In the ensuing discussion, therefore, we consider only the Hamiltonians -ffi.2 
and Hi,3 (Fig. ^), which describe interactions of "sites" in the x and y directions, respectively. 

For the wave function (|l9|) to be an exact eigenfunction of the Hamiltonian H, it is sufficient that the sixth-rank 
spinors 
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^Aipii^ipi ^l^^\2f^2^^2P2 (2)g^,x2 , (24) 

be eigenfunctions of the corresponding ceU Haniiltonians Hi 2 and -ffi,3. 

In general, when the site spinor vJ/^ai^^p jg j^q^ symmetric with respect to any indices, the possible states of two 
quartets of spins s = 1/2 consist of 70 multiplets. A wave function represented by a sixth-rank spinor contains only 
20 of them. Accordingly, the cell Hamiltonians Hi 2 and _ff 1 3 can be represented by the sum of projectors onto the 
50 missing multiplets: 

50 



Hia = J2^^Pk' 



fc=i 

50 



fc=i 

where the positive constants Afe and fik are the excitation energies of i?i,2 and Hi, 3, and the specific form of the 
projectors depends on 14 model parameters. 
Inasmuch as 

ffn,n+a|*s) = 0, i?„,„+b|*.) = 0, (27) 

for the total Hamiltonian (|2^) we have the expression 

H\^s)=0. (28) 

Consequently, is the ground-state wave function of the total Hamiltonian H, because it is a sum of nonnegative 
definite cell Hamiltonians. Also, it can be rigorously proved (see the Appendix) that the ground state of H is 
nondegenerate. 

As mentioned above, the specific form of the projectors depends on 14 model parameters, and in general the cell 
Hamiltonians (|2^), expressed in terms of scalar products of the type • s^, (s^ ■ Sj){sk ■ s;), etc., have an extremely 
cumbersome form. We therefore consider a few special cases. 

WhfiH the site spinor '^l^^^'^P is a symmetric fourth-rank spinor Q^^^'^p (corresponding to the two-dimensional AKLT 
modeltj), only the quintet component out of the six multiplets on each spin quartet is present in the wave function 
(|l|). The sixth-rank spinors ( |24| ) and (E5|) are symmetric with respect to two triplets of indices and, hence, contain 
four multiplets with S* 0, 1, 2, 3 formed from two quintets. Consequently, the cell Hamiltonian {H12 and H13 
coincide in this case) has the form 

66 

Hi,2^Y.^kPl^'- (29) 
fc=i 



If we set Afe = 1 {k — 1, 66), we can write Eq. ( p9[ ) in the form 

Hl,2 - ^4(81 + S2) + [1 - P2(Sl)P2(S2)] , (30) 

where Si is the total spin of the quartet of spins s = 1/2 on the ith site. Si ~ si{i) + 82(1) + S3(z) + S4(i), and P/(S) 
is the projector onto the state with spin S = I. 

If the four spins s = 1/2 at each site are replaced by a single spin 5 = 2 and if the wave function ( p^ ) is treated as 
a wave function describing a system of spins S = 2, the second term in the Hamiltonian (^ ) vanishes, and we 
arrive at the Hamiltonian of the two-dimensional AKLT model: 

Hi 2 = P4(Si + S2) = —Si • S2 + —(Si • S2? + —(Si • S2)^ + ^— (Si • S2)''. (31) 
^ ^ 28 40 180 2520 

Another interesting special case is encountered when the system decomposes into independent one-dimensional 
chains. This happens if the site spinor '^i^p'^p reduces to a product of two second-rank spinors, each describing two 
spins 1/2. For example, 

^^P''P{si, S2, S3, S4) = ifi^-'isu S3)v''fis2, S4). (32) 
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FIG. 4. Pattern of independent singlet pairs (double lines). 



In this case the Hamiltonians Hi 2 and -ffi.a contain interactions of four rather than eight spins 1/2 and have the 
form (0). 

The simplest case is when the site spinor '^■^^^'^P is a product of four first-rank spinors: 



(33) 



Now the system decomposes into independent singlet pairs (Fig. and the total Hamiltonian of the system has the 
form 



S ) • S7" 



(34) 



where and s,- are the spins forming the singlet pairs. 



4. SPIN CORRELATION FUNCTIONS IN THE GROUND STATE 



We now look at the problem of calculating the norm and the correlation function of the model described by the 
wave function (|l9|). The expression for the norm of the wave function G = (5's|^s) has the form 



G = n (vl/^»>"(n) M.^^P^-^P^(n)) 5.„A„+.gp„M„,.5.;A:^,^5p: 



a„ A = {1, 2, 3, 4}, 



(35) 



where Rapap isa4x4x4x4 matrix. 

According to the selection rules for the projection of the total spin S*^, only 70 of the 256 elements in the expression 

(^^^■'^ P{n) ^^P'^P {n)^ are nonvanishing. Consequently, the matrix R also contains at most 70 elements. If we 

regard the elements of R as Boltzmann vertex weights, the problem of calculating the norm reduces to the classical 
70-vertex model. 

Since the exact solution for the 70-vertex model is unknown, numerical methods must be used to calculate the norm 
and the expected values. 

To calculate the above-indicated expected values, we carry out Monte Carlo calculations on 20 x 20-site lattices. 
As mentioned, the ground-state wave function of the model depends on 14 parameters and, of course, cannot possibly 
be analyzed completely. We confine the numerical calculations to the case in which the spinor '^^p-'^p depends on one 
parameter a: 



COS a • Q 



(36) 
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FIG. 5. Dependence of the spin correlation function {s3(l)si(2)) on the parameter a. 



where a £ [ — 7r/2; 7r/2] , the spinor Q^^^^p is symmetric with respect to ah indices, and 

^Am^P = <^>(S1)^M(52)^''(S3)¥'''(S4). (37) 

In this case we have a one-parameter model with two well-known limiting cases. One corresponds to a = 7r/4, for 
which ^^>^^P = ^A^i/p^ ^j^g system decomposes into independent singlet pairs (Fig. the other limiting case 
corresponds to a = (our model reduces to the two-dimensional AKLT model in this case, the spins at each site 
forming a quintet). 

In the given model there are four spins s = 1/2 at each site, and the enumeration of each spin is determined by the 
order number of the lattice site to which it belongs and by its own number at this site. The spin correlation function 
therefore has the form 

/,,(r) = (s,(n)-s,(n-|-r)). (38) 

In determining the spin structure of the ground state, however, it is more practical to consider the more straight- 
forward quantity F{r): 

4 

F(r)= ^ (s.(n).s,(n + r)) = (S(n).S(n + r)). (39) 

ij = l 

The function -F(r) is left unchanged by a change of sign of a. This invariance is attributable to the fact that the 
spinor (^A^^p-'^p — Q^p-^p does not contain a quintet component, so that all the functions of this spinor are orthogonal 
to all functions of the symmetric spinor 



[Q^'P'^^'P' 

for all A, [1, V, p and A', /i', v' , p' 



[a^p^p - Q^P^P^^ ^ 0. (40) 



4 

In addition, since the total spin operator S at a site commutes with = ^ S; ■ Sj , we then have 



It follows from Eqs. (|35|), (40), and ( }4l| ) that sin a and cos a enter into the norm and into the expected value 
(^'|S(n) ■ S(n -f r)|^') only in even powers, so that -F(r) is invariant against a change of sign of a. We note, however, 
that only the total correlation function, and not (r), possesses symmetry under a change of sign of a. This assertion 
is evident, for example, in Fig. ^ which shows the dependence of /31(a) on a as an illustration. 

Figure ^ shows plots of F(r) for certain values of the parameter a. In every case it is found that the correlation 
function decays exponentially as r increases, differing from the one-dimensional model in that the preexponential factor 
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FIG. 7. Dependence of the correlation length on the parameter a. 



also depends on r. Figure |^ shows the dependence of the correlation length Vc on the parameter a. The correlation 
length is a maximum at the point a = (two-dimensional AKLT model), decreases as a increases, and at a = 7r/4, 
when the system decomposes into independent singlet pairs (Fig. ^), it is equal to zero. With a further increase in 
a the correlation length increases and attains a second maximum at a = Like the correlation function F{r), 

the function ri,{a) is symmetric with respect to a. It is evident from Fig. [7| that the parameter a has two rangess 
corresponding to states with different symmetries. In the range |a| < 7r/4 the correlation function F{r) exhibits 
antiferromagnetic behavior: 

F{r) cx (-!)'■- +'-«e-l'"l/''% (42) 

whereas the spins at one site are coupled ferromagnetically, (si(n) • Sj(n)) > 0. On the other hand, in the range 
7r/4 < |a| < 7r/2 the correlation function F{r) is always negative: 

F(r) oc -e-l'-l/'-'^ (43) 

and all the correlation functions at one site are also negative (Fig. ^). 

These ranges have two end points in common, a = ±7r/4, where Tc = 0. Whereas a = tt/A corresponds to the 
trivial partition of the system into independent singlet pairs, the case a = — 7r/4 is more interesting. 

In this case we have 

and the matrix (^vj/-^ p vJ/^m^^p^^ which enters into the equation for the norm ( |35| ) and the expected values, is 
transformed to 
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FIG. 8. Dependence of the spin correlation function at one site on the parameter a. 

— 2 (Q'^' ^' ''' A^^'^'''^ + l^j^' p' j^p-'^p'^ 
The symmetry of the spinor Q^^^f with respect to all the indices leads to the relation 



^qX'p'u'p QXfiup'^ _ ^j^X'p'v'p' QXpvp^j _ ^ 



\X' p u' p' 



j^Xpvp 



Equation (E^) therefore acquires the form 



l^^X'p'u'p' ^Xt.up\^ ^ l^^yp'^'p' ^Am^pJ^ ^ (5aA' 



"pp ■ 

^From the equation for the norm ( p5| ) we then have 



n 



A„ + aA' '5/'- IV-/,' , — 2 



2Af^ 



(45) 



(46) 



(47) 



(48) 



The latter equation has been derived with allowance for the relation 6^,yi g^xg^i \i — 5xx'- 

We now calculate the expected value (\I'|si(n) • Sj(n + r)!?*}. If sites n and n + r are not nearest neighbors, 
(\I'|si(n) • Sj{n + r)\'$) decomposes into the product of the expected values 



(*|si(n) • Sj(n + r)|*) = 2^^'^'-^ ( ^f^^VV 



s,(n) ^'^'^'^''(n) 
X (^'^^"'^'"'''"'P"' (n + r) Sj(n + r) ^'" p" {n + r) 

X SxX'Spp'S,y^'Spp>6x"X"'Sp"p"'S,y"iy'"6p"p'" =0. 



(49) 



Consequently, for a — —t: /A all the correlation functions at non- nearest neighbor sites are equal to zero. But if sites 
n and n + r are nearest neighbors, the corresponding correlation function assumes the form 



(^'|s,(l) • Sj(2)\^) = 22A'^' ^^-^^^^'(1) s,(l) ^^'^''''(n) 
X ^*-^"V"'''"'p"'(2) sj{2) ■^^"P"''"P" {2) 
X 9vX" gu'X'" 6xx' Spp' 5ppi Sp"p"' (5^' V" ^p"p' 



(50) 
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The exact calculation of the latter expression yields the following results (Fig. ^) : 

z = l, 2, 4, j = 2, 3, 4, 





■s.(2)) 




■si(2)) 


(S3(l)' 


■si(2)) 



27 
256' 



It follows from Eqs. (pT]) that 



^ (s,(l).s,(2)) = -J, 

as in the case of independent singlets {a — 7r/4). It can also be shown that all the correlations functions at one site 
are equal to zero. 

to write the cell Hamiltonian Hi^2 in explicit form for a = — 7r/4, we introduce the notation 

11 = Si(l) +S2(1) +S4(1), rSi=S3(l), 

12 =S2(2)+S3(2)+S4(2), 1s2-Si(2), 

(52) 

hi = li • Si + I2 • S2, 
h2 = ll ■ S2 + I2 ■ Si. 



Accordingly, choosing Afc = 1 (fc = 1, 50), we can write the cell Hamiltonian H12 in Eq. (]2^) in the form 

Hl^2 = Pl/2{h)Pl/2{h)Pl{si + S2) + F3/2(ll)^'3/2(l2)/l3 

+ P3/2{h)Pl/2{h)hi + Pl/2{h)P3/2{h)h5, (53) 



where 



207 49 3 , , 1 , x„ , X 15, 1 , 

256 + 64^^ ■ + 64^^ ' + 16^^^ ' ' 32^ 



f 7 11, , (54) 

^4 = 7 - -rSl • S2 + -ll • S2 + - (ll • Si)(li • S2) + H.C. , 

i 1 T T 

h^^ - - -SI-S2 + -h • Si + 7 [(I2 • S2)(l2 • Si) + H.C.] . 

4 o 4 4 

The cell Hamiltonian iJi_3 has the same form (|5^ ) but with a change of notation according to Fig. ^: 

' ll ^ Si(l) + S2(l) + S3(l), / Si S4(l), 



I2 =Si(3)+S3(3)+S4(3), \S2=S2(3). 

Of special interest is the case corresponding to a = ±Tr/2. Unfortunately, exact expressions for the correlation 
function cannot be obtained in this case, but the Hamiltonian can be written in explicit form. Since the site spinor 
vj/^Mi^p does not contain a quintet component for a — ±7r/2, the wave function of two nearest neighbor sites ( p^ and 
( p5| ) will lack a component with S* = 3. A more detailed analysis shows that 19 multiplets are present in the wave 
function of two nearest neighbor sites. In this case, therefore, the cell Hamiltonian has the general form 

51 

H^,2^Y.^kPl^- (56) 

fc=i 

For a definite choice of \k in Eq. (|5^) the cell Hamiltonian assumes the form 

Hi^2 - ^^2(11 + Sl) + P2(l2 + S2) + Pi/2(ll)Pl/2(l2)Pl(si + S2) 

+ ^3/2(ll)A/2(l2)/'l/2(ll + Sl + S2) + P3/2(l2)A/2(ll)A/2(l2 + §2 + Si), (57) 
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where the notations ( |52| ) and (55) are used for nearest neighbor sites along the horizontal and along the vertical, 
respectively. 

Our results suggest that the spin correlation functions decay exponentially with a correlation length ^ 1 for an 
arbitrary parameter a. We also assume that the decay of the correlation function is of the exponential type for the 
14-parameter model as well, i.e., for any choice of site spinor vj/^^'^''. This assumption is supported in special cases: 
1) the partition of the system into one-dimensional chains with exactly known exponentially decaying correlation 
functions; 2) the two-dimensional AKLT model, for which the exponential character of the decay of the correlation 
function has been rigorously proved.Lj Further evidence of the stated assumption lies in the numerical results obtained 
for various values of the parameter in the one-parameter model. 



5. GENERALIZATION OF THE MODEL TO OTHER TYPES OF LATTICES 

The wave function (^, (|l^) can be generalized to any type of lattice. The general principle of wave function 
construction for a system of spins 1/2 entails the following: 

1) Each bond on a given lattice has associated with it two indices running through the values 1 and 2, one at each 
end of the bond. 

2) Each bond has associated with it a metric spinor with the indices of the ends of this bond. 

3) Each site of the lattice (a site being interpreted here, of course, in the same sense as in Sec. ||) with m outgoing 
bonds has associated with it an mth-rank spinor with the indices of the bonds adjacent to the site. 

4) The wave function is the product of all spinors at sites of the lattice and all metric spinors. 

It is obvious that each index in the formulated wave function is encountered twice, so that the wave function is 
scalar and, hence, singlet. 

The wave function so constructed describes a system in which each lattice site contains as many spins s = 1/2 as 
the number of bonds emanating from it. 

To completely define the wave function, it is necessary to determine the specific form of all site spinors. The 
coefiicients that determine their form are then parameters of the model. 

The Hamiltonian of such a model is the sum of the cell Hamiltonians acting in the spin space of the subsystem 
formed by the spins at two mutually coupled sites: 

H = Y^m,. (58) 

(io) 

Each cell Hamiltonian is the sum of the projectors with arbitrary positive coefficients onto all multiplets possible 
in the corresponding two-site subsystem except those present in the constructed wave function: 

i/,,, = ^ AfeP^'^ (59) 

k 

Then H.i^j\^s) = and, accordingly, s) = 0. 

Consequently, is an exact ground-state wave function. 

We note that any two lattice sites can be joined by two, three, or more bonds, because this does not contradict the 
principle of construction of the wave function. Moreover, the general principle of construction of the wave function is 
valid not only for translationally symmetric lattices, but for any graph in general. As an example, let us consider the 
system shown in Fig. |^. The wave function of this system has the form 

VP, = (l)^^^^^^-^P^ (2)vI/P^-^-i (3)*^^-^ {^)9X^X.9^.^,.9.^u,9p,p,9r,r, (60) 

and describes a system containing ten spins 1/2. 

If the given lattice has dangling bonds (as occurs for systems with open boundary conditions), the resulting wave 
function represents a spinor of rank equal to the number of loose ends. The ground state of this kind of system is 
therefore 2'-fold degenerate (where I is the number of loose ends). For an open one-dimensional chain, for example, 
the ground state corresponds to four functions — one singlet and three triplet components. For higher-dimensional 
lattices this degeneracy depends on the size of the lattice and increases exponentially as its boundaries grow. 
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FIG. 9. Example of a graph corresponding to the wave function (|60[). 



6. CONCLUSION 



We have proposed a method for the construction of an exact wave function for a class of two-dimensional spin models. 
In general this model depends on 14 parameters, and its Hamiltonian is written as the sum of the Hamiltonians of 
nearest neighbor spin quartets. The exact ground-state wave function of the total system is also the exact wave 
function of each cell Hamiltonian. Since 20 of the 70 multiplets of two nearest neighbor quartets are present in the 
exact wave function, the cell Hamiltonians are the sums of the projectors with positive coefficients onto the other 
50 multiplets. These coefficients are the excitation energies of the corresponding multiplets. Different values of the 
coefficients correspond to different Hamiltonians. In this case, however, the ground-state wave function itself and the 
spin correlation functions in the ground state are identical for all Hamiltonians. This means that the ground-state 
wave function, as defined by Shastry and Sutherland,L3 is superstable. 

We have carried out Monte Carlo calculations of the spin correlation functions in the ground state for the special 
case of a model that depends on one parameter. For all values of the parameter the spin correlation functions 
decay exponentially with distance despite the complicated dependence of the correlation functions of nearest neighbor 
spins on the model parameter. It is justifiable to expect the spin correlations to decay exponentially in the general 
14-parameter model as well. 

In closing, we are pleased to acknowledge Prof. M. Ya. Ovchinnikova for helpful discussions of the problems treated 
in the article. This work has received financial support from the Russian Foundation for Basic Research (Grants 
No. 96-03-32186 and No. 97-03-33727) and from the Program for Support of Leading Scientific Schools (Grant No. 
96-15-97492). 



APPENDIX 

In Sec. I we have constructed the singlet wave function 

= n*'"'^"''"''"(n)3-.A„+.5p„p„,. (A.l) 

n 

for a system of 4M^ spins s = 1/2 on a square lattice. The following Hamiltonian was specially chosen for the resulting 
wave function: 

H = J2 ^"."+- + J2 Hn.n+h, (A.2) 



for which the wave function (A.l) is the zero-energy ground state: 

H\^s) = 0. (A.3) 



We now show that the ground state of the system is nondegenerate, i.e., the wave function satisfying Eq. ( [A.3| ) is 
unique. 
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In asmu ch as the Hamiltonian ( [A.2| ) is a sum of nonnegative definite cell Hamiltonians, any function satisfying 
Eq. ( [A.3| ) must satisfy all the cell equations 



(A.4) 



This means that Eqs. (A_^) and ( A.4 ) are equivalent 



We prove the nondegeneracy of the ground state of the Hamiltonian (A. 2) as follows. We first write the general form 
of the wave function for t he sy stem in question. We then determine the general form of the wave function satisfying 
one of the cell equations (A.4). Making note of the conditions imposed on the general form of the wave function 
by each cell equation and, at the same time, simultaneously satisfying these conditions for all the cell equations, 
we obtain th e gen eral form of the wave function satisfying all the equations (A.4) and, hence, satisfying the total 
Hamiltonian (A. 2). 

Any wave function of the given system can be written in the form 



(A.5) 



where the summation is over the 4Af^ indices A^, /i^, Vi, pi, c{\p,i'p) denotes coefficients that depend on these indices, 
and (^^i^^i'^iP3 (j) are arbitrary fourth- rank site spinors (in general, spinors at different sites can differ). 
We require that the wave function (A^) obey the cell equation 



^fn,„+a|*) = 0. 



(A.6) 



By the construction of the singlet wave function 




, which is matched by the cell Hamiltonian H^.n+a, any wave 



function at sites n and n + a that satisfies condition (|A.6|) is a linear combination of the 64 functions contained in 
the expression 



(A.7) 



because the cell Hamiltonian -ffn,n+a by definition is the sum of the projectors onto all multiplets [vp^nMni^nPn ('^^ 
and \l>-'>n+a/^n+a'^n+"iPn+a(n + a) arc definite site spinors occurring in the wave function (A.l)]. We note that these 64 



functions can be linearly dependent (as is the case, for example, for the two-dimensional AKLT model). 
Thus, the general form of the wave function satisfying Eq. (A.(;) can be written 



* = ^ c(AMz^p|z.„A„+a)g.„A„+.*^"^"""''"(n) 

XV]/An + E,A'n + E,l^n + aP„ + a(-jj _|_ J~J <J) l^j Pj ( j ) ^ 

j7^n,n+a 



(A.5 



where c(A/ii^p|fc'nAn-(-a) are coefficients that depend on the indices A^, pi, Vi, pi exclusive of the indices and An+a, 
and (^^i^^i^iPi i^]) are arbitrar y sit e spin ors. 

Comparing the functions (A.5) and ( |A.^ ), we deduce the followin g co nditions that must be met by the function 
(A.5) to obtain the general form of the wave function satisfying Eq. (A.6): 



1. The spinors at sites n and n -I- a must coincide with the site spinors of the wave function (A.l) 



^A„+aMn + al^„ + aPn + a,(-jj + a) = V]/^>» + '»P>» + ^'''" + aPn + a _|_ 

2. The coefficients c{\pvp) have the form 

C{\pvp) = c(A^i^p|i^„An+a)3i/„A„+a- 

From the equation 



I*) = 



we deduce analogous conditions on the general form of the wave function (A.5): 

$A„Mn--nP„(n) ^ ^A„Mn!^„Pn(n)^ 
<JjA„+bPn + b!^n + bPn + b^j^ + b) = Vj/ -^n + b /^n + b '^n + b Pn + b ^ 

C{\p.vp) = c{Xpi^p\pnPn+h)gp„p„+i,- 



(A.9) 

(A.IO) 
(A.ll) 

(A.12) 
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The simultaneous satisfaction of all the cell equations (A.4) requires consolidation of the conditions imposed by 
these equations on the general form of the wave function ( A. 5 ). Combining these conditions in succession, in each 
step we obtain the general form of a wave function satisfying the equations corresponding to these conditions. Upon 
satisfying all the condition s, w e obtain the general form of the wave function satisfying all the cell equations (A.4) 
and, hence, satisfying Eq. (A. 3): 



■9. - 



J2 c{Xniyp\Xniyp)]]_^' 



(A.13) 



where c{\pLV p\\pLV p) = c is a consta nt. 

Comparing the wave functions (A.l) and ( A.13|), we readily perceive that, to within an arbitrary factor, the 
general form of the wave function satisfying Eq. ( |A.3| ) coincide s wit h the wave function 'fg. Consequently, ^ s is the 



nondegenerate ground-state wave function of the Hamiltonian (A. 2) 
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